Random Walks, Electric Networks and 
The Transience Class problem of Sandpile 



Ayush Choure Sundar Vishwanathan 
Department of Computer Science and Engineering 
Indian Institute of Technology, Bombay 
ayush, sundar@cse.iitb.ac.in 

January 21, 2013 



Abstract 

The Abelian Sandpile Model is a discrete diffusion process defined on graphs 
(Dhar |15) . Dhar et al. [16j ) which serves as the standard model of self- organized 
criticality. The transience class of a sandpile is defined as the maximum number of 
particles that can be added without making the system recurrent ([!]). We develop 
the theory of discrete diffusions in contrast to continuous harmonic functions on 
graphs and establish connections between standard results in the study of random 
walks on graphs and sandpiles on graphs. Using this connection and building other 
necessary machinery we improve the main result of Babai and Gorodezky (SODA 
2007,(2]) of the bound on the transience class of an n x n grid, from 0(71^") to 
O(n^). Proving that the transience class is small validates the general notion that 
for most natural phenomenon, the time during which the system is transient is small. 
For degree bounded graphs, we demonstrate the first constant factor approximation 
algorithm for the transience class problem, based on harmonic functions. In addition, 
we use the machinery developed to prove a number of auxiliary results. We give 
general upper bounds on the transience class as a function of the number of edges 
to the sink. We exhibit an equivalence between two other tessellations of plane, the 
honeycomb and triangular lattices. 

Further, for planar sandpiles we derive an explicit algebraic expression which 
provably approximates the transience class of G to within 0{\E{G)\). This expres- 
sion is based on the spectrum of the Laplacian of the dual of the graph G. We 
also show a lower bound of Q,{n^) on the transience class on the grid improving the 
obvious bound of f2(n^). 



1 Introduction 

The abelian sandpile model (ASM) is a type of diffusion process defined on graphs which 
is closely related to the chip firing game investigated by Bjorner, Lovasz and Shor [7] 
and Tardos [33]. Indeed some of the results in the model, proved by Biggs [3], have 
an analog in the ASM of Dhar fTSl. The model proposed by Dhar has been studied 
in depth by the statistical physics community for investigating the phenomena known 
as self-organized criticality in the dynamics of sandpile formation. This formulation is 
known to be closely related to other interesting albeit diverse phenomena such as stress 
distribution in earthquakes, size distribution in raindrops, path length distributions in 
loop-erased random walks, for instance. For a nice overview, see the recent comprehensive 
survey article by Dhar [T3]. The ASM, though easy to define, has a very profound 
behavior and is far from being completely understood. Research in this area stretches 

'Extended draft appears in |12 | 
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across numerous disciplines such as probability theory, algorithmics, theory of computing, 
combinatorics, non-linear dynamics, fractals, cellular automata, to name a few. These 
connections have been beautifully summarized by Kleber in [^5] . Dhar [T3] also discusses 
some generalizations of the ASM like the Abelian Distributed Processors (ADP) model 
which is used to model a grid of abstract state machines along with many theoretical and 
practical applications. 

In the standard sandpile model, "sand particles" are added at the vertices of a (multi)graph. 
A site (vertex) is stable as long as the number of particles at the site remains less than its 
degree. Adding more particles would render the site unstable and is accompanied by the 
unstable site's passing a particle along each edge to its neighboring sites. This relaxation 
process is referred to as toppling. One of the sites known as the sink cannot topple. To 
ensure that every relaxation process eventually stabilizes, one needs the condition that 
the sink is reachable from every other site. As the system evolves, the sandpile goes 
through a a sequence of configurations. Those which can be revisited in any toppling 
sequence are called recurrent, the remaining ones are termed transient. Typically, one 
starts with the empty configuration and as particles are added, one moves through tran- 
sient configurations till a recurrent configuration is reached. Thereafter the configurations 
stay recurrent. The steady state behavior of a sandpile is characterized by its set of re- 
current states. It has been observed by physicists that for most natural phenomena, the 
time taken to reach a recurrent state is small. Hence any acceptable model must refiect 
this tendency to reach steady state rapidly and it becomes important to study the time 
taken to reach recurrence in these models. 

The essential parameter in our discussion is the number of particles which ensure 
recurrence. If particles are added randomly then a simple coupon collector type argument 
demonstrates polynomial bounds on the expected time to recurrence (as already mentioned 
in [2]). The other scenario is to add particles adversatively so as to avoid a recurrent state 
for as long as possible. This problem was highlighted by Babai and Toumpakari |4j where 
they define the requisite number of particles as the transience class of the sandpile. This 
later motivated the insightful work by Babai and Gorodezky [2] on grid based sandpile 
which are the most studied objects as compared to any other graph class because of their 
outstanding significance in statistical physics. In their path-breaking paper, Babai and 
Gorodezky [2J show that for the standard nxn square grid based sandpile, the maximum 
number of particles one can add before hitting a recurrent state is 0{n^^). This is a 
remarkable result in view of the fact that some closely related sandpile (for example line 
graph based) have transient state paths of length exponential in graph size. They use 
intricate combinatorial arguments based on particle conservation and the symmetry group 
of grid graphs to demonstrate the above mentioned bounds. However, simulations suggest 
a bound close to 0{n*) for the grid sandpile. Also the questions raised require analysis 
of the problem in a more general setting. 

Our contribution: We begin by showing a strong connection between the tran- 
sience class problem of sandpile and random walks on the underlying graph. Using LP 
duality and basic relaxation properties, we derive bounds on the transience class of a 
sandpile in terms of harmonic functions over the underlying graph. Similarly, bounds on 
sandpile impedances across any pairs of sites are obtained. These results form the core 
of our arguments which contrast the discrete model, sandpile, with the continuous ver- 
sion; random walks. For degree bound sandpile, we use the independent set properties of 
nodes with zero heights to demonstrate an algorithm which approximates the transience 
class up to constant factors. The algorithm works by computing harmonic functions over 
graphs and significantly tightens the connection between sandpile and random walks on 
the underlying graph. We then prove some basic properties of sandpile analogous to basic 
results in harmonic function theory, for example occurrence of worst case behavior at the 
boundary, reciprocity properties among any pair of sites, and so on. We derive and use a 
triangle inequality of potentials. This inequality provides sufficient fiexibility in analyzing 
the growth rates of harmonic functions at the cost of loosening the bounds. We use it 
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to obtain a bound on the corner to corner potential response on a grid network. Using 
some symmetry property of grids, we prove the main result of our paper which improves 
the bound on transience class from 0{n^'^) (by Babai and Gorodezky [1]) to 0{n'). We 
demonstrate a very general bound on the transience class in terms of sandpile size and 
the number of connections to the sink. We also show that in the case of planar sandpile, 
there exist explicit algebraic expressions which bound the transience class values. These 
are based on the spectrum of the Laplacian of planar dual of the given sandpile graph. We 
derive the expression for grid sandpile and leave it as a (somewhat technical) conjecture 
to establish bounds on its value. We believe that these would yield bounds as low as 
0{n'^) for this problem. In the last section we discuss some important and interesting 
open problems that would be of interest to the theory community. Our main contribution 
in this paper is to bridge the gap between discrete diffusions on graphs and the theory 
of harmonic functions on graphs. Indeed, random walks, electric networks, graph spectra 
and LP duality have been central tools in theoretical computer science. We hope that 
this paper initiates a theory of discrete diffusions analogous to the celebrated theory of 
mixing of Markov chains. 

1.1 Related Work 

Random walks and Sandpile: We begin by sketching a picture depicting an intuitive 
connection between sandpile models and random walks on graphs. On a graph G fix two 
vertices s and t. A simple random walk with a specified starting vertex v involves at each 
step, a choice of a neighboring vertex uniformly at random (See BoUobas [8] for a nice 
introduction). The potential sT^ti^) associated with v, and s and t as poles, is defined 
as the probability of reaching t before s starting from v. These tt functions, discussed 
at length in the next section are of paramount importance in analytic potential theory 
(see, for instance [37]). With site s as the designated sink, add particles at site v and 
observe the requisite number needed before a particle reaches site t. For any site which 
is ready to topple, if we label some particular particle among the set that are going to 
flow out, then the probability that it lands up at a particular neighbor, is uniform among 
the neighbors. If we add just enough particles at say Ny^ so that exactly one particle 
reaches t, then the probability of it being any particular particle is uniform and the path 
it takes from v to t looks just like the ones constituting s7rt(w). Informally speaking, 
any particle in the starting pile at v starts a random walk at v which terminates at t or 
s, whichever is encountered earlier. Intuitively, one expects that the probability s7rf(w) 
would be proportional to the reciprocal of iV„ with the proportionality factor accounting 
for discreteness and storage at sites, features absent in the usual network theory axioms. 
Our main theorem formalizes this connection and we derive as corollaries some properties 
of sandpile which are discrete analogues of the corresponding properties of random walks. 
Electric Networks: The classical theory of electric networks along with the well under- 
stood connections with random walks ([30], [17],[2I]) has some very powerful and intuitive 
results. These results have recently found applications in almost every important area of 
theoretical computer science. Christiano, Kelner, M^dry and Spielman [llj have recently 
announced the fastest known algorithm for computing approximate maximum s — t flows 
in capacitated undirected graphs. Using the electric current flows in this network with s 
and t as poles, their algorithm constructs approximate fiows. Earlier Kelner and M^dry 
|22| used arguments based on random walks to formulate the fastest known algorithm 
for generating spanning trees from uniform distribution. Spielman and Srivastava |33| 
construct good sparsifiers of weighted graphs via an efficient algorithm for computing ap- 
proximate effective resistance between any two vertices, a result which is quite insightful 
on its own. The list of important results which use harmonic functions in an essential 
manner goes on. The benefit of this confluence of research in different classical areas is 
indeed mutual. For example, in their path breaking paper, Arora, Rao and Vazirani [Tj 
give an 0{^/\ogn) approximation algorithm for computing graph conductance. Our goal 
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has been to use the theory of harmonic functions in analyzing sandpile behavior (in the 
context of diffusion) in analogy with the theory of random walks on graphs. The results 
we report in this paper open up the possibility of analyzing those properties of ASM 
which may not have been possible using purely combinatorial arguments. 
Other Results on Sandpile: As already mentioned, research problems on the abelian 
sandpile model span across numerous areas. Recent advances with a complexity theoretic 
flavor include proof of the one-dimensional sandpile prediction problem in LOGDCFL 
by Peter Bro Milterson [28J. Also, Schulz [31J mentions a related NP-complete problem. 
The group structure of the space of recurrent configurations, first introduced by Dhar, 
Ruelle, Sen and Verma in [TB], is also a fertile area of analysis. Cori and Rossin [13J show 
that sandpile groups of dual planar graphs are isomorphic. Toumpakari [38J discusses 
some interesting properties of sandpile groups of regular trees where questions related to 
group rank are studied and the paper is concluded with an interesting conjecture on the 
rank of all Sylow subgroups of the sandpile group. Specific families of graphs like square 
cycles C^, x C„, 3 x n twisted bracelets, etc have been analyzed. We refer the reader 
to [21], [12|, [20|. 

2 Preliminaries 

2.1 Introduction to the Abelian Sandpile Model 

Our notation and terminology follows Babai and Gorodezky j^. 

Definition 1. A graph G is an ordered pair {V{G)tE{G)) where V{G) is called the set 
of vertices and E{G) is a set of 2— subsets of V ^ possibly with repeated elements, the set 
of edges. 

This is referred to as a multi-graph in literature but we will use graph for brevity. The 
degree of a vertex v ^ V \s defined as the number of edges in E which contain v. Two 
vertices v and u are called adjacent (or neighboring) if {u, v) £ E. A path between two 
vertices u and v is an ordered sequence of edges ei, 62, . . . , such that u e ei, w e 
and for all values of i, H e^+i 7^ 0. The graph G is connected if there exists a path 
between any pair of vertices. 

To model an Abelian Sandpile Model, we take a connected graph G with a special 
vertex called the sink, denoted s € V. Non-sink vertices in G are called ordinary vertices 
and this subset will be denoted by Vo — V — {s} . 

Definition 2. The configuration of a sandpile G is a map c : Vo — N, which will be 
represented as a vector. The weight of c is |c| — J2vev ''(^)- 

The configuration c records the number of sand particles contained in each of the 
ordinary sites. The empty configuration is the zero vector. The capacity of a site is the 
maximum number of particles that it can hold and is one less then the degree of the node. 

Definition 3. An ordinary node v is said to be unstable in a configuration c if c{v) > 
degree(v). The configuration c is said to be unstable if any site under it is unstable, else 
it is referred to as stable. 

When a site is unstable it is said to topple, that is it passes on some of its particles to 
its neighbors. When a site v topples once, it loses degree(v) particles and each neighbor 
of V acquires a particle for every edge common with v. The sink node never topples. 
Starting with the empty configuration, we keep adding particles one by one on sites of 
our choice and topple them when necessary. 

The ASM evolves in time through two modes, particle addition at sites and relaxation 
of unstable sites via topplings. A toppling sequence is an ordered set of configurations 
where every configuration can be obtained from the previous one by toppling some unsta- 
ble site in it. Note that the event of many sites becoming unstable simultaneously poses 
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no complication since the order in which they are subsequently relaxed does not affect the 
final stable configuration that is obtained at the end of toppling sequence. Elementary 
proofs of such confluence properties can be found in the pioneering paper on ASMs by 
Dhar [15J. See also Babai and Toumpakari (4j. 

Notation: We write ci > C2 if Vw,ci(w) > C2 and ci h C2 if there is a toppling sequence 
which takes ci to C2. Finally we write, ci — > C2 if 3c3 > ci such that C3 h C2. We say that 
a configuration C2 is reachable from ci if ci C2 and unreachable otherwise. In other 
words, one can add particles to certain sites in ci so that there exists a toppling sequence 
leading to C2. Note that reachability is transitive, i.e. ci — >■ C2, C2 — >■ C3 ci — > C3. 

Theorem 2.1. (\l ^ .f ?j ) Given any configuration c, there exists a unique stable configu- 
ration a{c) such that c h a{c), independent of the chosen toppling sequence. 

Property 2.1. //c h cr(c), then kc h ka-{c) 

Associated with every toppling sequence is the count on the number of times each site 
has toppled, the vector of toppling potentials, also referred to as the score vector in [2]. 
These toppling potentials are very closely related to the electric potentials that develop 
at various nodes when power source-sink are appropriately applied, a connection which 
we will discuss in detail in the coming sections. 

Definition 4. Assuming ci h C2, the toppling potential function z^^'"^^ : Vo ^ N is 
defined as z'^^''^^(v) : the number of times v toppled in a toppling sequence from ci to 02- 
We denote z^^*^^^^ by z'^. 

This function is well defined as the number of times a particular site topples is inde- 
pendent of the toppling sequence chosen, already noted in |2j. The proof employs the fact 
that the principal minor of a connected graph's combinatorial laplacian is of full rank. 

A configuration is called recurrent if it is reachable from any configuration. As already 
mentioned, we say that a configuration Ci is reachable from a configuration Cj if by adding 
some particles to Cj (possibly at multiple sites) and subsequently relaxing it, we can obtain 
Ci. A configuration is transient if it is not recurrent. The set of recurrent configurations 
is therefore, closed under reachability. 

Property 2.2. If 3c' such that there is a toppling seqquence from c' to c in which every 
site has toppled at least once, then c is recurrent. 

The proof follows from the fact that the existence of such a toppling sequence precludes 
the existence of forbidden sub-configurations and hence makes the configuration recurrent. 
For a complete discussion on forbidden sub-configurations and recurrence of configurations 
the reader is referred to [15J and [14J. 

We analyze the process of adding one grain at a time to the sandpile and study its 
evolution. As in the standard theory of Markov chains, recurrence characterizes the 
long term (steady state) behavior of sandpiles. Our investigation is concerned with the 
maximum number of particles that can be added while staying transient. Following Babai 
and Gorodezky [2], for a sandpile S we define. 

Definition 5. The transience class of S denoted by tcl(5), is defined as the maximum 
number of particles that can be added to S before reaching a recurrent configuration. 

In view of propertv l2.21 we can bound the transience class from above by the maximum 
number of particles that can be added before all the nodes have toppled at least once. 
Showing that this bound is tight upto constant factors is also not very hard. We defer a 
fuller disussion of questions of this nature to our subsequent manuscript. 
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2.2 Basics of Harmonic Functions and Potential Theory 



For a very nice introduction to harmonic functions on graphs, we refer the reader to the 
beautifully written paper by Benjamini and Lovasz [5J and to Teles [37J for a thorough 
view. We start with some important definitions and fundamental properties. Given a 
connected graph G and a function tt : V{G) — ^ R, we say that tt is harmonic over Vh if, 



The remaining vertices (lying in V" — Vh) are called the "poles" of tt. The set Vh is 
also called the interior of tt with vertices adjacent to the set of poles referred to as the 
boundary. We see that the value of vr at any vertex in Vh is the average of its value in 
the immediate neighborhood. In case of multi graphs, we take the appropriate weighted 
means, where the weights are the number of common edges. This leads us to the first 
basic property, 

Property 2.3. Any non-constant harmonic function can assume its extreme values only 
at the set of poles. 

It follows that every non-constant harmonic function has at least two poles, its maxima 
and minima. Such functions are completely determined by their values on these vertices. 
Formally speaking. 

Property 2.4. Uniqueness: If two functions harmonic on Vh agree on the boundary, 
they agree everywhere in the interior. 

More generally, we have the following property, 

Property 2.5. Given a set of poles, a harmonic function is uniquely determined modulo 
scaling and translation by a constant. 

Properties 12.41 and 12 . 51 important as they allow one considerable freedom in construct- 
ing harmonic completions of functions defined over the boundary set. This problem is 
the discrete analogue of the classical boundary value problems in complex analysis. We 
will describe two important examples in which these function arise naturally. 

Random Walks on Graphs: Consider a graph G and two special vertices s and t. The 
potential associated with u, with s and t as poles, s'^ti'v) is defined as the probability 
of reaching t before s starting from v. One can check that the function tt so defined 
is indeed harmonic on the set V — {s,t}, with the maximum value of 1 at the node t 
and the minimum value at s. The generalization to the multi-pole situation is also 
straightforward . 

Electric Networks: Consider a resistive electric network (i.e. a circuit made up entirely 
of resistors). Let gT^tiv) be the potential that appears at node v when unit potential is 
applied across t and s. Using the equation of charge conservation (KirchofF's node law), 
one can show that these potentials are harmonic on all nodes except s and t. 

The main implication here is that one can intuitively think of the electric network the- 
ory as an analysis of random walks of electrons on the underlying graphs. Consequently, 
results from network theory can be used to prove interesting facts in other related ar- 
eas. As an example, consider the problem of constructing the harmonic completion of a 
function with given boundary values. All one needs to do is to take the corresponding 
circuit and apply potentials equal to the boundary values on the boundary points. The 
potentials that will appear on other nodes can be computed using basic linear algebra (the 
only non-trivial step involves inverting the combinatorial Laplacian of G) thus allowing 
construction of harmonic completions efficiently. We outline below three very basic and 
fundamental results of network analysis which will be needed in the following sections. 




(1) 
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Theorem 2.2. Superposition Principle: The superposition principle states that for 
all linear systems, the net response at a given place and time caused by two or more stimuli 
is the sum of the responses which would have been caused by each stimulus individually. 

Theorem 2.3. Compensation Theorem: If the impedance Z of a branch in a network 
in which a current I flows is changed by a finite amount dZ , then the change in the 
currents in all other branches of the network may be calculated by inserting a voltage 
source of —IdZ into that branch with all other voltage sources replaced by their internal 
impedances. 

Theorem 2.4. Reciprocity Theorem: In its simplest form, the reciprocity theorem 
states that if an emf E in one branch of a reciprocal network produces a current I in 
another, then if the emf E is moved from the first to the second branch, it will cause the 

same current in the first branch, inhere the emf has been replaced by a short circuit. Any 
network composed of linear, bilateral elements (such as R, L and C) is reciprocal. 

The reciprocity theorem can be restated in terms of just potential sources and poten- 
tial measurements using the notion of effective resistances between pairs of nodes. The 

effective resistance between a pair of nodes u and v, R^j f{u. v) is defined as the potential 
difference which develops between u and w if a unit current source is applied across u and 
v. 

Lemma 2.1. Potential Reciprocity Lemma : If taking s and t as poles with 7r(s) = and 
7r(f) = 1 induces a potential o/s7rt(w) at node v and interchanging the roles of v and t 
induces sT^v {t) at t then, 



Proof : Consider the given network G with the special node s G V{G). We refer 

to the corresponding modified network G{e) obtained from G by adding an edge with 
resistance 1/e between every node and s. In particular, G(0) = G. Furthermore, we refer 
to an edge between s and u by s~u. We will be using the current source version of the 
reciprocity theorem. If applying a unit current source across st results in a potential of 
V across 6~v, then applying unit current source across sv results in a potential of v units 
across st. The value of this potential v can be expressed, using Ohm's law, as the ratio 
of current through the edge and the resistance of the e— edge between the particular node 
and sink. Since both potentials are equal in magnitude, we can say that on G(e), 



This follows from observing that applying a unit current source across sv is equivalent 
to applying a voltage source of Reff{s, v) across sv. Because of linearity, it follows that a 
potential of sTT-u{t)Reff{s,v). appears at node t. Similarly so for the other configuration. 

This equation holds for arbitrarily small values of e. Consequently it holds for graph 



In particular, when the effective resistances across s and t are the same as s and v, 
we have sT^t{v) =« 7r„(f). In the following discussion, we will omit the left subscript (s) 
from sTTt whenever it is clear from context. We say that a walk P is an instance of sTt* 
if it starts at some vertex v, avoids s and ends at t. The following lemma may already 
be known to experts. Since we could not find it in literature, we present it with a simple 
proof. 

Lemma 2.2. A triangle inequality for potentials 



Refs{s,t)s-Kt{v) = Reff{s,v)s1Ty{t) 



(2) 



Reff{s,t)s'Kt{v) = Reff{s,v)sny{t) 



G(0). ■ 



■Ki{j).Trj{k) < TTi{k) 



(3) 
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Proof : Let 'P{k,j) be the set of all walks from k to j avoiding s. This set can 
be partitioned into two components, namely the walks passing through i and the ones 
avoiding it, denoted by Vi{k,j) and Vi{k,j) respectively. For any walk P which is an 
instance of tt^, let the probability of occurrence be ■Ki{P). Then by definition, 

^Ak)= E ^^■(^)+ E ^^■(^) 

Similarly, 

Using these two relations, we obtain 

Mj).n,ik)= E ^^p)- E '^^■(^)+ E ^^(^)- E ^^■(^) (4) 

P£V{j,i} P&Vi{k,j) PeV(j,i) P&Viik,]) 

Consider the first term on right side in equation Being a probability measure, the 
value of X]pe-p(j i) '^ii-P) bounded above by 1. Every s-avoiding walk from k to j passing 
through i can be decomposed into two components, a walk from k to i avoiding j and a 
walk from i to j. This implies, J2peVi{k,j) '^AP) = Y.p^VJ{k,^) ^^iP)-J2pev{t.j) ^^P) < 
^pgp^^^ TTii^P). The first term therefore has the following bounds, 

E n^ip)■ E ^^(^)^ E ^»(^) (5) 

Pev{j,i) Pev,(k,j) Pev-j{k,i) 

For bounding the second term, observe that any s-avoiding walk from fc to i which 
passes through j, can be treated as a juxtaposition of a walk from k to j, avoiding i, and 
a walk from j to i. Hence, 

E MP)- E ^^(^) = E ""^p) (6) 

PePU,i) PeVi{k,j) PeVj{k,i) 

Using equations |4l [5] and [6] we get 

^,:(j)-7r,(fc) < E ^'^(^)+ E MP) = Mk) 

PeVjik,i) PeVj(k,i) 



Remark: The utility of this inequality becomes clear when interpreted in the context of 
electric networks. Consider a network such that the node with ground potential is fixed 
and we are allowed to apply power at any other node and observe the resulting potentials. 
The inequality implies that if applying a potential Vi at i produces unit potential at node 
j and applying V2 at node j produces unit potential at node fc, then applying V1.V2 units 
at i produces at least unit potential at node k. 

3 Reducing the transience class problem to estimating 
harmonic functions over graphs 

We first consider the single site particle addition strategies. We will later show that the 
effect of allowing particle addition at multiple sites on our transience class estimates is 
inconsequential as far as our estimates are concerned. 
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Definition 6. Consider a sandpile S with nodes u and w. The sandpile impedance of 
the ordered pair {v, w), Rs{v, w) is defined as the maximum number of particles that one 
can add at v before a toppling at w occurs. 

Note that unlike the impedance of electric networks, sandpile impedance is not sym- 
metric in its arguments, i.e. in general Rs{v,w) ^ Ra{w,v). To estimate its value, we 
introduce the following LP relaxation. 

max 

< ^ z{v') - d{v).z{v) < d{v) - 1 

v' 

\/u^v:0<'^z{u')-d{u).z{u) < - 1 

z{w) < 0,z > 0,x > 

The values of Xy (the number of particles added at v) and z (the vector of toppling 
counts) that are realized above are a feasible solution of this LP and hence the optimum 
of this LP yields an upper bound on the Rs{v, w). With the fixed sink node. ,s. we define 
77^ (t;) as the potential at node v when a unit potential is applied at node w. In terms of 
these potential functions, the following bound holds. 

Lemma 3.1. The optimum value of the above LP is bounded from above by the following 
value, 

-J^Y,(diu)-l).7rUu) (7) 

Proof: We consider the following relaxed version of the given LP. 

max Xv 

z{v') - d{v).z{v) + Xy < d{v) - 1 

yu^v: z{u') - d{u).z{u) < d{u) - 1 

ziw) < 0,z > 0,a; > 

From the weak duality for LPs, it follows that to obtain an upper bound of a on the 
optimum value of the above system, it suffices to find a feasible solution of the dual LP 
of value a. The dual is the following : 

min ^(d(M) - l).Y{u) 

U 

^ Y{u') + Y' - d{w).Y{w) > 

u' r^W 

yuy^w: Y{u') - d{u).Y{u) > 

Y{v)>l,Y>0,Y' >0 
Consider the following set of equations 

J2 Y{u') + Y' - d{w).Y{w) = (8) 

yuy^w: Y ^("') - d{u).Y{u) = (9) 

Y{v) = 1 
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A non-negative set of values satisfying the above set is feasible for the dual LP. We 
find these by considering the resistive circuit S, obtained by replacing each edge in S by 
a unit resistance. We assign ground potential to the sink, and inject current at node w 
such that it gets unit potential. The potential that develops on any node u is s^^wiu)- 
The potential value at node v, 'Kw{v), can be used to scale the input current at w thereby 
scaling all the potentials as well, such that potential at node v becomes unit. It follows 
that the values Y{u) = TTwiu)/nyj{v) and Y' equaling the value of the current injected 
form a feasible solution of the dual LP. The objective value at this point is, 

^{d{u) - l).TTyj{u) 



This yields an upper bound on Xy . To obtain a lower bound, consider the complemen- 
tary problem of finding x'^, the minimum number of particles that must be added at v to 
observe a toppling at w. The following LP's objective value forms a lower bound on x'^, 

minja;^} 

< ^ z{v') - d{v).z{v) -h < < d{v) - 1 

Vu ^ u : < ^ z(m') - d(u).z{u) < d{u) - 1 

z{w) > l,z > 0,a; > 
The proof of the following lemma is analogous to the previous case. 

Lemma 3.2. The optimum value of the above LP is bounded from below by the following 
value, 

(10) 



Proof: Consider the relaxed version of the above LP, 

min{<} 
J2 z{v')-d{v).z{v)+x'^ > 

Mu^v: ^ z{u') - d(u).z{u) > 

z{w) >l,z>0,x>0 

From the weak duality for LPs, it follows that to obtain a lower bound on the optimum 
value of the above system, it suffices to find a feasible solution of the dual LP. The dual 
is the following : 

max Y{w) 
"^Y{u')+Y' -d{w).Y{w) < 

yuy^w: ^K) - d{u).Y{u) < 
Y{v) <1,Y>0,Y' >0 
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Consider the following set of equations 



Y{u') + Y' - d{w).Y{w) = (11) 

yu^w: r(u') - d{u).Y{u) = (12) 

Y{v) = 1 

As before, any non-negative set of values satisfying the above system is feasible for 
the dual LP, and therefore forms a lower boudn on the objective value. We find these 
by considering the resistive circuit S, obtained by replacing each edge in by a unit 
resistance. We assign ground potential to the sink, and inject current at node w such 
that it gets unit potential. The potential that develops on any node u is sT^w{u). The 
potential value at node w, 7rtu(w), can be used to scale the input current at w thereby 
scaling all the potentials as well, such that potential at node v becomes unit. It follows 
that the values Y{u) = T^wiu) / it w(v) and Y' equaling the value of the current injected 
form a feasible solution of the dual LP. The objective value at this point is tTw{v)^^. ■ 

Clearly the maximum number of particles that can be added at v before toppling some 
w is just one less then the minimum number that need to be added at v to topple w, 
that is — Xv + I- Using equations and ([TU]), the following two-sided bounds are 
obtained. 



-^-l<.„<-^^(d(.)-l)..4-) (13) 

Further, define the potential profile of the circuit when unit potential is applied at 
node w as 

rs(u;) =^(d(?;)-l).^„(i') 

V 

Using this notation, Rs{v, w) = Xy satisfies the following general bounds. 

Lemma 3.3. Rs{v,w) is 0{rg{w).Tryj{v)^^). 

Lemma 3.4. Rs{v,w) is n{TT^{v)~'^) . 

To find the maximum number of particles one can add at v before every other site 
topples, one simply needs to consider the maximum value of Rs{v^w) over all values of 
w € Vo- Consequently, one can find the maximum number of particles that can be added 
at a single site before every other site topples, by considering the maximum of Rs(v,w) 
over all pairs (Vjw). This value, iiia.'Ky^yj{Ts{'w)Tr^{v)~^}, also forms a bound on tcl{S) 
as allowing particle addition at multiple sites above gives the same estimates. This fact 
follows from essentially the same line of argument that was used for finding the upper 
bounds except that in this case instead of Xy, the objective function to maximize is Xu 
where Xu is the number of particles added at u. 

Theorem 3.1. tcl{S) is 0{ma.Xy^yj{ryj{S)TTyj{v)~^}). 

As each of the tt^{v) lies between and 1, the value of r^{S) is therefore bounded 
between 1 and |2i<^(5')|. Hence, we have the following relaxed upper bound on the value 
of tcl{S), 

Lemma 3.5. tcl{S) is 0{\E{S)\.ma-Ky^y,{'Kw{v)~'^}). 
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These results quantify the relationship between sandpiles and random walks on graphs. 
However, the theory of potential functions on graphs boasts of several very intuitive and 
beautiful results, e.g. the recipocity theorem. In a later section, we show that the 
parallelism between sandpiles and electric networks runs deeper by demonstrating the 
sandpile versions of some well known basic results in network theory. 

We start by considering some simple properties of potential funtions. For example, 
consider the propertv 12.31 which says that the maximum and minimum of potential func- 
tions occur at poles. An elementary proof by contradiction is easily conceivable. In the 
case of sandpiles, one can think of a similar notion of maxima/minima in terms of ease of 
percolation of particles. When we add particles anywhere (may be more than one site), 
then is it so that the last site to topple will be adjacent to sink? Or consider the dual 
problem. We are allowed to add particles at one site only. For any particular site whose 
toppling we wish to delay for as long as possible, is it true that the best strategy is to 
add particles at a site adjacent to sink? To rephrase, is it true that for any site w, the 
value of Rsiv, w) is maximized for some v adjacent to the sink? 

Note that both these questions are two sides of the same coin in case of random walks 
because of the reciprocity properties discussed above. The first has a direct analogy for 
sandpiles. 

Lemma 3.6. The last site to topple is always adjacent to sink. 

Proof : Assume that the last site to topple is not adjacent to sink. Since each of its 
neighbors has already toppled, it has received at least as many particles as its degree and 
has become unstable at least once contradicting the assumption that the particular site 
has never toppled. ■ 

We now fix the site under observation and ask the same question about the site where 
we add particles. 

Lemma 3.7. For a given w, the estimate of Rs{v,w) is maximum when v is at boundary. 

Proof :For a fixed w, the value of Ts{w) is fixed. One has to show that the value of 
n.u}{v) is minimum for some vertex v adjacent to sink, s. This clearly follows from the fact 
that for every internal node u, t:{u) is a convex combination, in particular the weighted 
arithmetic mean, of the 7r(.) values at its neighbors. This means that t:{u) is bounded 
between the values spanned by the neighbors, so it cannot be an extreme point. ■ 
Note: This lemma talks about the estimate and not the exact value of Rs{v,w). The 
lemmas prove that while using the Theorem 13. 1[ it is enough to consider both sites on 
the boundary set (i.e. adjacent to sink). The following lemma is the sandpile analogue of 
the classical potential reciprocity lemma from network theory. 

Lemma 3.8. Sandpile Reciprocity lemma : If adding p particles at v causes toppling at 
w then adding 2\E{S) \ -P particles at w causes a toppling at v. 

Proof : Using theorem 13. 3( the ratio of Rs{v,w) to Rs{w,v) can be bounded. 

Rs{w,v) ^ maxRsjw.v) ^ Tr^jw)-^ .r{v) ^ 2\E(S)\ ^"ff'^'"'^'^ 
Rs{v,w) ~ mm Rs{v,w) 7r„(w)-i ~ 'Reff{w,s) 

Where the last inequality follows from the potential reciprocity mentioned in Lemma 
12.11 Given p — Rs{v, w), we get the required bound on Rs{w, v) in terms of p. 
U 

One can go even further by using the fact that the maximum value of Reff is |^(s)| 
(attained for paths) and minimum value is at least 1/|_B(S')| (attained for just a pair 
of adjacent nodes with many parallel edes betwene them). Babai and Gorodezky [2] 
conjectured the following. 
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Conjecture 1. (J2I) Assume that for sandpile X: the induced sub-graph on the set of 
ordinary vertices is connected. Then the transience class of x ( the largest weight of any 
transient configuration ) is the height of the tallest transient stack of grains placed on a 
single site. 

The conjecture is equivalent to saying that using single site particle addition strategies 
one can attain the transience class bounds. However this is not so. A counter-example 
by Sunic appears in [3]. We present a simpler counterexample and an intuitive reason 
why this conjecture is false in general. Assume that in a sandpile x with finite transience 
class, one can attain the bounds by adding particles at the single site v and w topples last. 
After the last particle is added, adding one more is supposed to make the configuration 
recurrent. Which means this heaviest transient configuration, when relaxed should have 
every site filled to its maximum capacity, except for w which has not yet toppled and 
contains particles less then the maximum capacity. If there exists some other site which 
is not filled up to maximum capacity, one can add particles there and fill it up (only till 
it stays stable, of course). Hence the validity of conjecture rests on the rather unlikely 
premise that in the heaviest transient configuration, every site but one is filled to its 
maximum stable capacity. A condition which one would think unlikely when there are 
no symmetries in x (trivial automorphism groups). We will present an example with 
non-trivial symmetries to demonstrate that even in this case, one cannot expect such a 
strong property. 

Consider the grid sandpile Xn with n = 4. We add particles at the top left corner. As 
expected, the last site to topple is the bottom right corner. Here is the stable configuration 
corresponding to the heaviest transient configuration with a stack of particles placed on 
the top left. 



All the sites which have fewer then 3 particles, can be topped up without inducing a 
toppling at the bottom right corner node, and then adding a particle at top left corner 
induces recurrence. This demonstrates that single site particle addition do not work for 
this pair. We discuss a related open question in the section of future work. 

4 A constant factor approximation for transience classes 
of degree bounded graphs 

We will now argue that in the case of degree bound graphs, using some combinatorial 
properties of independent sets in sandpile graphs along with arguments similar to the 
ones outlined in the previous section for the lower bound derived in Lemma 13.21 one can 
derive tight (up to constant factors) lower bounds for the transience classes. 

In the sandpile S, let x'^ be the minimum number of particles that can need to be 
added at node ti such that the node w topples at least once. Clearly x'^ = Xy + 1. As in 
the previous section, the following LP relaxation forms a lower bound on the value of x'^ . 



3 3 3 



3 3 2 



3 3 2 3 



2 3 2 



mm X, 



V 




degree(w) — 1 



Vu 7^ u : < 2, •^('^O ~ degree(w).z(M) < degree(u) — 1 



z{w) > l,z > 0,x^ > 
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Consider the scenaxio in which we add x'^ particles at v and allow the configuration 
to settle down to stability. Furthermore, let the number of particles that appear at any 
node u in the resulting stable configuration be The following lemma gives a lower 
bound on the particle count x'^ in terms of these height functions 

Lemma 4.1. The value of the particle count x'^ is bounded from below by the following 
value, 

-l^EM-)-.^ (14) 

Proof: Clearly < /i„ < degree(w). In the linear program stated above, replacing 
each pair of constraints of the type < J^v'^^v ^i"^') ~ d{v)-z{v) + Xy < d{v) — 1 by 
'I2v'~v •^('^') ~ d(v).z{v) + Xy = hy, maintains the property that the optimum value is at 
least a lower bound to the exact solution to the particle count x[,. The altered LP is the 
following, 

min x'^ 

z{v') — degTee{v).z{v) + x'^, = h{v) 

\/u V : z{u') — degree(u).z(u) = h{u) 

z{w) > 1,-2 > 0,4 > 

To bound the optimum value, we will consider the dual of this minimization program 
and find a suitable feasible point. The value of the cost function at that point will be 
used as bound. The dual is the following maximization program. 

max Y' + Y^h{u).Y{u) 

U 

^ F(u') + i"'-degree(w).y(w) < 

u' 

Mut^w: ^ Y{u') - degree(u).F(u) < 

Y{v) < i,r > o,r' > 

As in the previous lemma, we consider the following set of equations, whose feasible 
region is inside the one corresponding to the dual we mentioned above. 

Y{u') +Y' - d{w).Y{w) = 

\/u w : ^(^') — degTee{u) .Y (u) = 

Y{v) = 1 

A non-negative set of values satisfying the above set is feasible for the dual LP. We 
find these by considering the resistive circuit S corresponding to the graph S (with each 
edge having unit resistance) . The sink node is assigned ground potential and just enough 
current is injected at node w so that it attains unit potential. In terms of the standard 
potential functions described earlier, the potential that develops on any node u is sT^w{u). 
The potential value at node v, Trw{v), can be used to scale the input current at w thereby 
scaling all the potentials as well, such that potential at node v becomes unit. It follows 
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that the values Y(u) = tTii;{u)/t:w{v) and Y' equahng the value of the current injected 
form a feasible solution of the dual LP. The objective value at this point is, 

^ ' u 

Since the current, y , is a positive quantity, the above derivation implies the Lemma. □ 
We will now show that the lower bound is at most a constant factor smaller then the 
upper bound. We start with defining an auxiliary set of variables hy such that, 



hy = 1 — hy > 1 

From the positivity of 7r(.) function and the domination relation hy < hy, the following 
inequality follows. 

r = ^ E ^«-»(") ^ ^ E '^"-»(") (15) 

We will now bound the value F from below. 

Lemma 4.2. (Dhar \14^ , Babai and Gorodezky J^) In any stable recurrent configuration, 
for every edge, both the incident vertices cannot have zero particles. 

This follows from the fact that in any stable recurrent configuration, the last toppling 
of one of them would have taken place after the other and so the second node necessarily 
has at least one particle. For every edge u, v, at least one of hu and hy is > 1. 

Corollary 4.1. The set of nodes I = {v : hy = 0} form an independent set. 

Assume further that the graph satisfies (A), i.e. the maximum degree is A. Using 
this we will obtain the required bounds. Consider any vertex v £ Vu along with its 
neighborhood (see figure ([Ij). The function 7r(.) is harmonic over this neighborhood, so 
we have 

degree(u)7r(w) ~ ^ 7r(M) 
from which we get the sum of 7r(.) over any neighborhood in V/,. as. 



■n{N{v)) ~ tt{v) + tt{u) — (degree(u) + l)7r(?;) 

The contribution of any local region in the restricted potential trace F depends on 
just two possibilities regarding hy. We deal with them separately. 

- Case hy = 0: v £ I, consequently none of its neighbors are in /. In this case the 
contribution to F is the sum E«^i,7r(w) = degree(w)7r(u) = dtlfll{v)+i '^ ^'")) ■ 

- Case hy = \: v ^ I, and some of its neighbors are in /. Assume the worst case 
scenario when all the neighbors are in /. The contribution to F in such a situation 

is just 7:{v) = dcgree(.) + l ^(^(«))- 
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Center node selected Center node's neighbors selected 



Figure 1: Neighborhood around vertex v; two cases, v selected and not selected in / 

It follows that any neighborhood N{v) contributes at least ^^^^.^^^^■^^j^ 7r(-/V(f )) to F, 
regardless of the specific values of hy. Since the degree is bounded by A. this translates 
to a minimum contribution of -^^7r{N{v)). Therefore, the value of F in equation 
for any independent set /. is at least 



(A + l)7r^(«) 

The upper bound yielded by Lemma [3. II is bounded from above by, 



Y^TT^U) (17) 



Using equations (fT7|) and (|T6| . we obtain the following two sided bounds for degree 
bounded graphs, 

1 (A - 1) ■c-^ 

,Tryj{u)<Xy< _ ^^.^ } ^TTi^{u) (18) 



(A + l)7r^(u)^ 7r^(w) 

Theorem 4.1. For any sandpile with bounded vertex degrees, the minimum number of 
particles that need to be added at any vertex v to observe a toppling at any vertex w is 
equal, up to constant factors, to the following expression, 

^ ^n^{u) (19) 



TTwiv) 



The computation of transience class tcl(S') requires evaluating the above expression 
for all possible combinations oi v,w S V{S). Computing the function 7r(.) can be done 
in very efficiently following the recent path breaking work by [23], [Ml) [34] on solving 
symmetric, diagonally-dominant linear systems. Consequently, finding the pair with worst 
estimates is also easy to do. 

Corollary 4.2. There exists a polynomial time algorithm which computes the transience 
class of a degree bound sandpile up to constant factors. 



5 The case of Grid Sandpile 

As noted in the introduction, the sandpile associated with the n x n grid is of particular 
importance. We define it formally below. 

Definition 7. Consider the nxn grid graph. Attach an extra sink node to the boundary 
such that there is a single edge to each non-corner boundary node and double edges to 
the corner nodes. We denote both the sandpile and the corresponding circuit by GRID„. 
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Notation: For the purposes of labeling the nodes, we assume the grid is embedded 
canonically in the first quadrant of with a corner coinciding with (1, 1). Every node 
on the grid is labeled with the coordinates it occupies in the lattice. The labels are 
(*5 j)i 1 < J < f^- The sink node is labeled s. 

Babai and Gorodezky [2J have shown that td(GRID„) — 0(n^^ '^°^^). In this section we 
will improve this bound to O(n^). The following is a broad outline of our proof of Theorem 
15.11 We bound the potential profile r(GRID„) and min^j^tu t^v{w) separately and estimate 
the bound on td(GRID„) using Theorem 13.11 The bounds on r(GRID„) are obtained 
using ideas based on charge conservation, along the lines of the classical Ampere's Law of 
electrodynamics. For bounding the value of min^^^ n^lw), we show that values of 7ri,(c) 
and 7rtu(c) (where c is the center) can be used to obtain estimates on 7r„(z«). Using grid 
symmetries we prove monotonicity properties which imply that the minimum value of 
7r^(c) is obtained when w is a corner. Finally we bound the value of 7ry{c) by constructing 
a harmonic function with power applied at corner such that unit potential appears at the 
center. The construction of this distribution uses a certain potential domination property 
of the center over edges and the fact that the grid graph can be expressed as the Cartesian 
product of paths. The amenability of paths in constructing harmonic distributions and 
the classical superposition theorem (Theorem l2.2l) play key roles in the construction. We 
begin in the next subsection, with the potential domination property. 



5.1 A potential domination property of the center 

We will consider the case when potential is applied at a corner and prove a kind of potential 
dominating property of the center over the corner opposite to the power source. The 
proofs of these monotonicity properties require concepts involved in proving convergence 
properties of iterative algorithms which solve boundary value problems. This procedure 
is known as the Jacobi Method 0. One starts with assigning the given values to boundary 
points and zero to every other node. In every iteration, the value of any internal node is 
updated according to the values of neighboring nodes just after the preceding iteration 
ended. When the linear system is irreducible weakly diagonally dominated (as in our case), 
it produces a set of values converging to the final solution. For a proof of convergence 
we refer the reader to |35]- This technique is folklore in basic finite element analysis 
and belongs to the much more general class of algorithmic constructions of solutions to 
Poisson's equation. The speed at which the values converge to the solution is intimately 
tied up to the rate of mixing on the underlying graphs. We will show that when potential 
is applied at a corner, the values that appear on the nodes in any iteration obey a 
simple monotonicity property, thereby implying that the solution (which is the point of 
convergence of these points) obeys the same monotonicity property. 

Definition 8. Corner Monotonicity: Let / be a function defined on a finite n x n grid, 
/ : Z„ X Z„ — >■ N. We say / is corner monotone with respect to (1,1) if f{p) < f{q) for 
any pair of lattice points p and q such that the segment q — p is either perpendicular to 
the diagonal passing through (1,1) or along one of the edges passing through it and q is 
closer to the diagonal then p. 

Corner Monotonicity with respect to other corners is defined likewise , see figure ([2]). 

Let Co be the starting set of values with 1 assigned to (1,1) and to every other node. 
Let ct be the set of values resulting from iteration number t. ct+i is obtained from cj 
using the following conditions of harmonicity of functions. For any node u. 

Lemma 5.1. If Cq is corner monotone, then Ct is corner monotone for all values oft. 



^See the wikipedia entry for the Jacobi method 
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Figure 2: Corner Monotonicity with respect to (1, 1) 



Proof: Without loss of generahty, assume that q > p (the other possibihty wiU be 
iniphed by symmetry) . First consider the case when the segment p — g is perpendicular 
to the diagonal through (1,1). So if g is of form {x, y) then p is (a; — 1, y + 1). 

We will proceed by induction on the number of steps of the algorithm. Before the 
first iteration, time t = 0, corner monotonicity of cq is trivial. Assume C( is corner 
monotone, we need to show that ct+i{q) < Ct+i{p)- We now use equation ((20)) . Note 
that by induction hypothesis, each term in the expression of Ct+i{q) is dominated by the 
respective term of C(+i(p) which implies ct+iip) ^ ct+i{q)- A special case arises when p 
lies on the diagonal itself. Here we make use of symmetry of the grid. When p lies on 
the diagonal, its northern neighbor is mirror image of eastern neighbor and likewise for 
southern and western neighbors. The eastern and southern neighbors are common with 
q. The remaining two of g's neighbors are dominated by these two. Again, by induction 
hypothesis, the inequality follows. The remaining reasoning is same as the standard case. 

The other case of p — g being parallel to an edge through (1, 1) edge is analogous. ■ 

The limiting value of c is the harmonic distribution that results when a unit poten- 
tial is applied at the node (1,1). It satisfies the same monotonicity properties that the 
distributions c(t) satisfied, for all values of t. This gives us the following lemma. 

Lemma 5.2. When a potential is applied at a corner, then the resulting potential distri- 
bution is corner monotone with respect to that corner. 

Using this, the following potential domination property of the center can be inferred. 

Lemma 5.3. When potential is applied at a corner, the potential at the center of the grid 
is higher than at any site on the opposite boundary. 

Proof: Using corner monotonicity, we claim that when power is applied at node (1, 1) 
and unit potential is observed at some node {(n, i)} on the opposite edge , then the center 
of the grid (f , f), also has at least a unit potential. The reasoning behind this assertion 
is as follows. Because of symmetry, the site (i,n) also has at least unit potential. On 
the line connecting these two sites, say L, the potentials first increase till one reaches the 
intersection with the diagonal Di : x = y and then decrease monotonically. This follows 
from the corner monotonicity lemma as the starting configuration is corner monotone. So 



18 



both (1, 1) and L n Z?i(aiid in case LO Di is not a lattice point, the two points closest to 
it) have at least unit potential. Assume there exists a point on the line segment joining 
(1, 1) to LnDi, say k whose potential is less then unity. Then every point to its right has 
potential less then unity, following corner monotonicity. Similarly for every point right 
above it. But this two sets partition the circuit into disjoint pieces, one of which contains 
(1, 1) and other contains L D Di. Any random walk starting from L n Di and ending at 
(1, 1) has to pass through this set. The potential that appears on L H Di cannot exceed 
the maximum value taken by any point in this set. This contradicts the assumption that 
potential at i H Di is greater then that at k. Hence such a k cannot exist implying that 
every site on the line joining (1, 1) and L H Di has at least unit potential. ■ 
Note: Lemma 15.31 can be rephrased in the following manner. If applying a potential of 
p{n) at a corner produces unit potential anywhere on an opposite boundary node, then 
applying p{n) at any corner is enough to produce at least a unit potential at the center. 
We will later see an example of a harmonic distribution with a single positive pole at a 
corner and unit potential at some point on the opposite edge. The utility of this Lemma 
lies in the fact that in general constructing harmonic functions with an arbitrary pair 
of poles and known value at some arbitrary point is not easy. In our case, we need the 
potential that appears on a corner when potential is applied at the opposite corner. Our 
efforts so far, to construct a distribution with a pole at corner and known response at the 
opposite corner, have been fruitless. However, using Lemma 15.31 in conjunction with the 
triangle inequality for potentials, we obtain fairly good estimates of the corner to corner 
potential correlations. We believe that the estimates we obtain are close to the square of 
the true value. 



5.2 The case of corner to corner 

We will now obtain a lower bound on the minimum value of TTyiw) for any pair v and w. 
We will show that minimum values of 7r„(w) are obtained when both v and w are points 
on the boundary of grid. Let the center of the grid be denoted by c. Then, using Lemma 
12.21 (triangle inequality of potentials), we obtain 7r^,(c)7rc(u) < ■k^{w). Using Lemma |2.1|, 



7^^„(c) = ■Kciw).—^ 

Reff{s,w) 

Clearly, /3 min^, 7rc(u)^ < min„.„ 7r„(z/;), where j3 is the minimum, upto constant factors. 



value of over all possibilities of boundary nodes w. The following lemma bounds 



the value of /3. 

Lemma 5.4. The value of /3 defined above is lower bounded by some constant. 



Proof: The bound is derived in two parts. We first derive a lower bound on the 
numerator. Consider any node w. The effective resistance between sink node s, and 
w decreases if we reduce any edge's resistance. This follows simply from the Reiligh's 
monotonicity principle. We reduce all the resistances, except the ones incident on w, to 
zero. This effectively leaves only node w connected with s by 4 parallel edges. The net 
resistance of this configuration is 1/4. This is an absolute lower bound on the effective 
resistance between any node and sink (as the argument is independent of the location of 
node). 

An upper bound on the value of denominator follows from the fact that it is a parallel 
combination of a unit resistance with a network. The net resistance of parallel combina- 
tion of ri and 2 is at most min{ri, r2}. 

The above two facts give the required bounds on the value of /3. ■ 
To show that TTciv) is minimum when w is a corner node, we need another potential 
monotonicity lemma, the sandpile analogue of which appears in Babai and Gorodezky [2]. 
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However, the use of monotonicity properties in our proof is essentially different. Babai and 
Gorodezky [2J use monotonicity along with the pigeonhole principle based combinatorial 
arguments to derive bounds on tc?(GRID„). These arguments are first made on the infinite 
grid. Using monotonicity, [2] bound the region which particles touch when they are added 
to single sites and by ensuring that the sizes of the regions are small, one can assume 
that the boundary is not touched and pretend to be on the infinite grid itself. Our use of 
monotonicity is much more straightforward in the sense that we want to find the pair of 
vertices with the worst estimates and monotonicity properties lead us directly to them. 

Definition 9. Center Monotonicity: Let / be a function defined on a finite n x n grid, 
/ : Z„ X Z„ — > N. We say / is center monotone if /(p) < f{q) for any pair of lattice points 
p and q such that the segment q — p is aligned perpendicular to some axis of symmetry 
and q is closer to it then p. 

Lemma 5.5. When potential is applied at the center, then the resulting potential distri- 
bution is center monotone. 

The proof of above lemma is completely analogous to the previous case. The center 
monotonicity lemma implies that if we apply a potential at the center, then the corner sites 
have the lowest potential (among all non-sink nodes) . Rephrasing in terms of reciprocals 
of 7r(.), we get the following upper bound on the maximum value of 7r„(z«)^^ over all pairs 
w and V. 

Lemma 5.6. // applying the potential p{n) on a corner induces unit potential at the 
center, applying K .p{nY / j3 at any node induces unit potential at every non-sink node, 
where K is a constant. 

The only remaining information is the value, p{n) , of the potential that when applied 
at a corner, induces a unit potential at the center. In the next section we will see an 
example of a such a harmonic distribution. 

Note: In the preceding discussion, we have conveniently assumed that n is odd, else no 
such center site would exist. It is however easy to extend the discussion to the case of 
even n. 

5.3 Constructing a harmonic distribution over GRID„: determin- 
ing the corner to center response 

Our current goal is to construct a harmonic distribution with power applied at a corner 
such that at least unit potential appears at the center, or in other words bound the 
value TTu^^^^^^ (c). In general, constructing harmonic distributions with arbitrary poles 
and known values at some node in general is difficult. However, in our present problem, 
we will bound this quantity using the fact that the grid is the Cartesian product of paths 
and that potential functions on paths are easy to construct. 

Consider the (path) line circuit which has n nodes. The last node is connected to 
ground potential through a unit resistance. We apply a potential of n + 1 units through 
a unit resistor at node labeled n and observe that unit potential appears at the corner 
vertex (labeled 1). See figure 



3 2 1 

Figure 3: Line Circuit 

A harmonic distribution on GRIDn'. Now take the n x n grid with ground connection 
attached to each of its boundary nodes on the left and bottom edges through unit resis- 
tances. Power sources are applied at the top and right edges through unit resistances. 



n + l 
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We apply a potential of (n + l).i at the boundary nodes (n, i) and {i,n). At the special 
corner node (1, 1), we apply + n. One can check that the potential that appears at any 
grid node is V{i,j) = i.j. In particular, unit potential appears at node (1,1), i.e. 
V{1,1) = 1. This construction is a particular case of constructing a harmonic distribution 
on the Cartesian product of two graphs given a harmonic distribution on each of them. 
The generalization is discussed in the full version. 

Using the superposition principle (Theorem 12. 2p . the potential value at (1, 1) due to 
these 2n — 1 power sources is the sum of potential values that would have appeared when 
these power sources would have been used one at a time with all other sources short 
circuited, at their respective positions. Also, among all the nodes on the top and right 
edges, there exists one with the maximum potential response at (1,1), i.e. where when 
unit potential is applied, the potential at (1,1) is maximum. Again using superposition 
principle, if all the power sources are applied at this site alone with all other sites con- 
nected to sink, at least unit potential appears at (1,1). The value of this new power 
source is + n + 2.(" + 1).* = n^^ + = 0{n^). However, the site on which 

power source is applied has exactly one connection less with sink compared to the circuit 
GRID„ in which we apply power through unit resistances. To remedy this, we add an extra 
edge to the sink. One can check that this modification is non-essential and one needs to 
change the power applied at this corner by the amount of current flowing this new edge, 
which amounts to an at most constant factor change in the power applied. The same 
potentials appear at other nodes. Using reciprocity and (degree) regularity of ordinary 
sites in GRID„ one can interchange this power node on an edge and the corner node (1,1) 
to obtain the following lemma. 

Lemma 5.7. In GRIDn, applying 0{n^) potential at (1, 1) induces 0(1) potential at some 
point on the top edge. 

Using Lemmas 15.71 and 15.31 'we get the following result: 

Lemma 5.8. In GRIDn, applying 0{n^) potential at (1, 1) induces at least 0{1) potential 
at the center. 

Remark (Improving the lower bounds on tcl(GRIDn))'. Lemma 15.71 observes a pair of a 
vertices, both at boundary, such that applying 0{n^) potential at one of them induces 
at least unit potential at the other. The paragraph preceding this Lemma outlines the 
proof of this property by shifting the power sources to the best response point on the 
edge. If however one shifts these sources to the worst response point, the existence of 
a complementary vertex, with respect to the corner, can be proved. This pair has the 
property that applying 0{n^) potential at one of them induces at most unit potential 
at the other. Using the Lemma |3.4[ we obtain that the number of particles that can be 
added at one of them without toppling the second one is lower bounded by fl{n^). This 
is an improvement over the obvious lower bounds of fl{n'^). 

Corollary 5.1. tcl{GRIDn) = il{n^). 

5.4 Bounding the potential profile of GRID„ 

Using the discussion preceding Lemma [3.51 we can get a bound of O(n^) for r(GRID„), 
which yields an 0{n^ log n) bound on teZ(GRID„). Here we improve the bound on r(GRID„) 
to 0{n) using current conservation arguments and the regularity of normal nodes of 
GRID„. Using property (|3.6[) . one knows that the last site to topple is always at the 
boundary. Hence, when using theorem 13.11 for estimating transience classes, we know 
that the current source will always be added to a boundary node (adjacent to the sink). 
Consider our network, GRID„, with a current source attached to some node v adjacent to 
the sink, such that the potential of v is unit. Note that the total current flowing in, ?, is 
bounded above by the degree of v (4 in this case). To see this, consider the equation of 
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current conservation at v. v is at unit potential and each of the neighbors' potential is 
non-negative. Consequently, at most a unit of current flows through each incident edge. 
So the total outflow is bounded from above by 4 units. Hence, the total inflow from the 
current source is also bounded from above by 4 units. Thus, applying a current source 
of 0(1) units produces a unit potential at the node v. Using this fact, we can prove the 
following bound on r^(GRID„). 

Lemma 5.9. Consider any vertex v on the boundary of grid in GRIDn- The potential 
profile Ty{GRIDn) induced due to 0{1) current source at v is 0{n). 

Proof: Denote the potential at node v by ■n{v). The total current going into the sink 
node is, say, i and is equal to 

i = Tr{v) 

(D,s)e_B(GRID„) 

The current going in to the sink is a sum above all the constituting currents through 
each of the incident edges. Since the potential of s is zero, each of these currents is equal 
in magnitude to the potential of the neighboring nodes. If the set of nodes on corners are 
denoted by C„ , those on the interior of the boundary edges by /„ , and the union of these 
two by Bn, we can rewrite the value of i in the following form. 

z = 2. ^ n{v) + J2 ^i^) (21) 

■ueC, vein 

Since all potentials are positive, the following inequality follows from equation (|2ip. 

z > '^(") (22) 

Now consider the sequence of smaller (n — 2k) x (n — 2k) concentric grids nested in 
the larger n x n grid. Define the sets C„-2/c, In-2k and Bn~.2k analogously for each of 
these. For any element v G -^n-2(fc+i)j denote by n{v) the unique neighbor lying in In-2k, 
and for a v £ C'„_2(fe+i), denote the two neighbors by ni{v) and n2{v). Then, for each of 
these smaller grids, the net current entering through the set Bn-2k is zero. In terms of 
potential functions, the condition can be stated as 

Y (2.7r(w) -7r(ni(t;)) -7r(n2(w))) + ^ (7r(w) - 7r(n(w))) = 

'U£C„_2(fc + l) wG-fn-2(fe + l) 

Separating the vertices belonging to boundaries of different grids, we obtain 

2 E E ^(^) = E ^(^) (23) 

Again using the fact that all potentials are positive, we get 

Y Hv) < E ^(^) (24) 

i'e-B„_2(fc+i) ■ue-B„_2fe 

Every vertex belongs to the boundary of exactly one concentric grid. Using equation 
^ and dH]), we get 
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(25) 
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I) e GRID, 



lies, 



Since the degree of normal vertices is 4, we have the following bound on the potential 
profile of grid when power is applied at some vertex of boundary B„ ■ 



5.5 Transience Class of GRID„: a new bound 

Using Lemmas 15. 6[ 15.81 and 15.41 we obtain the following result which bounds the value of 
maXy,iu{Try{w))~^ from above. 

Lemma 5.10. In GRID„, applying 0{rJ) potential at any site induces at least Oil) 
potential everywhere. 

Using Lemma [5 . 1 i 1 5 . 9 1 and Theorem 13. II we have the following bounds on td(GRID„). 
Theorem 5.1. tcl{GRIDn) = 0{n'^). 

Remark: While the bounds proved above mark a substantial improvement over the current 
known 0(n'^® *'°^''^), experiments suggest a bound of somewhere O(n^). The estimates on 
the value of potential profile has little scope of improving substantially. Constructing 
the harmonic distribution with more care seems to be a plausible approach. Another 
possibility lies in exploiting the planarity of the sandpile graph. We will explore this 
avenue in further detail in the next section and obtain closed form expressions on bounds 
of tcl(S), when S is planar, in terms of the spectrum of the Laplacian of the dual of S. 

6 The case of planar Sandpile 

After showing the intimate relationship between the transience class of a sandpile and the 
harmonic functions over the underlying graphs, we will now show that if the underlying 
graphs are planar, the bounds on transience class can be expressed in a much more explicit 
algebraic form. Consider the sandpile S and the corresponding circuit, both of which will 
be assumed to be planar. 

In the circuit 5, we apply a unit potential across some boundary edge and observe 
the potential at some boundary node (as has already been noted, boundary nodes suffice 
for our worst case analysis). Now take the dual planar circuit of S, say S (for a detailed 
discussion of dualising operations in context of harmonic functions, see Benjamini and 
Lovasz |5j). For every edge in the original graph, there exists exactly one edge in the dual 
graph. Call these edges dual of each other. There is a special edge in the circuit, the 
power edge, across which the potential source is attached. Its dual edge becomes the unit 
current source in the dual circuit. The potentials at nodes in original circuit satisfied the 
Kirchoff 's current law (the condition of harmonicity of voltages) . The potential difference 
across each edge becomes the current flowing through the respective dual edge in the 
dual graph. And equations of Kirchoff's current law become those of loop law in the 
dual. Since all the currents satisfy the loop law, the potentials thus developed satisfy the 
current law as well. The estimation of potential difference across any edge in the first 
circuit is equivalent to estimating the current across the dual edge in the dual circuit. 
For any boundary vertex, its potential difference with that of sink equals in magnitude 




(26) 
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the current through the boundary edge incident on this node. The same current flows 
through the dual edge in the dual circuit. So, estimating the potential of a boundary 
node in the original graph is equivalent to estimating some current in the dual graph. 

As an example of interest, consider the grid graph with sink attached to border. We 
take its dual graph. See figure (|4]). The original circuit is shown in black lines and the 
dual in red lines. The dotted edges on the top right corner are the power sources of the 
two circuits. If for a unit potential applied at top right corner produces x at the bottom 
left corner (labeled a'), then unit current source through the top right edge produces a 
current x through the bottom left resistor (labeled ab). 



Potential Source - ■ 
Current Source . . 
Zero potential — 





planar circuit duality => Vgi - Igi^ 



Figure 4: Grid Circuit with its dual 

Note: The current source is placed across the dual edge. Which means that if the dual 
edge is connecting the nodes u and v, such that source is attached to u and sink to v, then 
a unit current fiows from v to u in the edge uv, to maintain flow conservation equations. 
In network theory terms, the edge uv itself is the current source, and as such its internal 
flows must not be taken into account while writing the Kirchoff 's equations. Consequently 
we delete the edge uv and simply attach a current source at u and sink at node v. The 
graph obtained after deleting edge uv, dual to edge e in 5, is called the restricted dual 
of S and denoted by Se- We will use the same symbol to denote the underlying graph. 
The current source is attached to vertex u and sink to vertex v. Let the combinatorial 
Laplacian of this graph be denoted by L, the potentials that appear at each of the nodes 
because of the current flowing by vector Z and let / be the vector containing net currents 
flowing in at any node. The equations of Kirchoff 's laws at each node can be succinctly 
written as. 



LZ ^ I 

The vector I has all entries except for a 1 at position corresponding to node u and 
— 1 for node v. Given the values of L and /, we need to estimate the potential difference 
between nodes p and q, or equivalently, the current in the edge pq. Because matrix L 
is singular, it is not possible to resolve the question by the usual methods of estimating 
certain entries of the inverse matrix. However, using the fact that L is symmetric one 
can indeed almost invert it enough to suffice for our purpose. See, for instance the paper 
by Wu |39j to compute two point resistances in networks. 
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Lemma 6.1. Consider a resistive network S with Laplacian L whose eigenvalues are 
Ao = < Ai < ■ ■ • A„_2 < A„_i and ^' is the unitary matrix containing the eigenvectors. 
The i^^ column, tpi, is the eigenvector corresponding to Xi. If unit current is injected at 
node u and taken out from node v then the magnitude of current in edge pq is given by 



~ ' ^ Xk 

0<fc<n-l 



(27) 



Proof : Denote by L[e) the matrix L + eJ. Note that for e > 0, L(e) is invertible, 
unhke L. Call its inverse G{e). Denote the row of G{e) corresponding to node a by 
G(e)(a). 

Z{a) = lim(G(e)(a)./) 
Knowing the value / explicitly, we can write 

Z{a) = lim(G(e)(a, u) - G(e)(a, v)) (28) 

where G{e){x,y) is the entry in the row of the node x and column of node y. The 
formula for G(e) is, 

G(e) = *A(e)-i*t 

where ^ is the unitary matrix containing the eigenvectors of L (and consequently 
of L{e) and G(e)) as its columns, 'i'^ is its hermitian and A(e) is the diagonal matrix 
containing the eigenvalues of L{e). Exactly one of these eigenvalues is e (for connected 
graphs). The corresponding eigenvector has every entry 1/n. Let the eigenvalues be 
Aq = e < Ai < . . . A„_2 < A„_i. We obtain the following expression for G{e){x, y). 



Gie)ix,y)^^+ 5: t^M^M. (29) 
We ^-^ At + e 

where 4'k(x) is the entry of node x in the A:*'* eigenvector. Using (|28p and (1^9)1 we get, 



Z{a)^ y Ma)iMu)-Mv)V (30) 

Xk 

The amplitude of the current in the edge e between nodes p and q is equal in magnitude 
to Z{p) - Z{q). That is. 



E 



(V'fe(p) - ipk{q)){il'k{u) - ipkiv)y 



Xk 

0<k<n~l 



111 the circuit 5, the possible locations of a power source are just the set nodes con- 
nected to the sink. The possible current source and sink nodes in the dual S are exactly 
the ones corresponding to dual (boundary) edges in S. Let {e^lfLg be the set of boundary 
edges sandpile S. Denote by kipq-, the current in edge pq in Se^. when current source is 
applied across the edge dual to Ck in S. Using this notation with Lemmas 13.51 and 16.11 
we can bound the transience class of S purely in terms of eigenvalues and eigenvectors of 
the Laplacians of restricted duals of S. 
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Lemma 6.2. For planar sandpile S, tcl{S) = |S'|.0(maxfe maxpq |fc*pq^|)- 

Proof: The quantity minpg \kipq\ gives the smallest boundary node potential value 
when the power edge in the original graph S* is e^. Further, miufc miup^ \kipq\ gives the 
minimum boundary node potential value over all possible placements of power source. 
The reciprocal of this value is equal to the quantity maxi,^tu{7ri„(w)~^}. The potential 
profile of any circuit is bounded by the size of the graph. ■ 

Remark: Note that one could as easily invert the toppling matrix of the sandpile to work 
out a similar formula. Being the principal minor of a connected graph's Laplacian, it is 
invertible. But as opposed to Laplacians, their principal minors are not well explored. 
While we know the eigenfunctions of Laplacians of most common classes of graphs, the 
same cannot be said about the principal minors. It is in this light that the above sim- 
plification becomes important. Its power can be displayed by considering the example of 
the grid. We discuss this example further in the last section on future work. 

7 Sandpile with k connections to the sink 

We will illustrate the power of our reduction by proving the following theorem for any 
sandpile S having at most k connections to sink. 

Theorem 7.1. td{S) is 0{{\S\ + 2)'=) 

In a sandpile S, let the vertices connected to the sink be labeled vi,V2, ■ ■ ■ ,Vk. We 
have already seen that the worst case estimates of m.a.Xy^ui{'^wiv)~^} are obtained, when 
both the nodes labeled v and w are on the boundary, i.e. when both are directly con- 
nected to the sink. Henceforth, we will limit our discussion to only those pairs of sites 
which are from this boundary set Bs = {filfLi- So, tcl{S) is bound from above by 
0(|£;(5)|.maxi<,,,<fc{^„^(i;,)-i}). 

Consider the circuit S. The sink site, s, has k edges incident on the remaining graph 
as shown in figure([5]). Assume that among all ordered pairs of vertices {v,w) from Bs, 
the value of n^iv) is minimum for the pair {vi,Vj). We will obtain lower bounds on the 
value of TTy. (vj) in terms of size of S and fc, the number of connections with the sink. 




Figure 5: Replacing the circuit with a reduced circuit 

We have already mentioned that 7r„ . {vj ) is the potential that appears at node Vj when 
unit potential is applied at node Vi. This is equivalent to saying that the probability of 
a random walk starting at node Vj hitting node Vi before the sink is equal to TTy.{vj). 
Consider the graph S — {s} obtained from S by deleting s (this is the induced sub-graph 
over the set Vo). In this graph, if there is no path from node Vi to Vj then n^- (vj) is zero, 
consequently its transience class is infinite. This refers to the condition that the sandpile 
graph must stay connected even if we remove the sink, for the transience class problem 
to make sense, as noted by Babai and Gorodezky [2J. If every path between some pair 
of nodes passes through the sink, then one can add an unbounded number of particles at 
one of these sites without being able to induce a toppling at the other site. The following 
lemma summarizes this observation. 
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Our goal is to construct a sequence of graphs {S'p}™ with Sq = S and having some 
nice monotonicity properties on the value of TTy^(vj). In Sq — {s}, consider any edge, e, 
whose deletion does not disconnect Vi and Vj. We apply a unit potential at site Vi and 
observe the potential at node Vj. For the edge e, either reducing its resistance (eventual 
contraction) or increasing its resistance (eventual deletion) decreases tt^. (wj). This fact 
follows from a trivial extension of exercise //.4.15 in Bollobas [8J. 5*1 is obtained from 5*0 
by deleting/contracting e, whichever operation reduces tt^,. (vj). For the sake of uniformity, 
the node labels of 5*1 are inherited from 5*0 canonically. When an edge is contracted, the 
resulting node can be labeled with the label of any of the colluding nodes. Note that at 
no time nodes Vi and Vj get merged as this would increase tt^,. (vj) to 1, which contradicts 
the property that in any iteration, its value can not increase. 

Similarly, in the p*'' iteration, Sp^i is transformed to Sp by picking any edge of Sp-i — 
{s} whose deletion does not disconnect Vi^Vj and contracting/deleting as appropriate. 
Observe that in every iteration, the total number of edges goes down by one. So the 
algorithm terminates with some graph Sm which satisfies the following property. 

Property 7.1. If Sm is the graph obtained after the last iteration, deleting any edge in 
Srn ~ {s} disconnects Vi and Vj. 

Which is equivalent to saying that the graph Sm — {s} is just the path from Vi to Vj, 
depicted by the circuit on right in figure ([5]). Note that none of the edges incident to 
the sink get deleted/contracted during the whole process. Naturally, the path length is 
bound from above by the size of E{S) — k. 

Now consider the circuit Sm- We will be observe only the adjacent nodes of s. The 
resistance connecting any two of these (which occur consecutively on the path Sm — {s}) 
is bounded above by the total path length, |i?(S')| — k. Unit potential is applied at node 
Vi and the node s is maintained at zero potential. The embedding shown in figure(l5]) 
demonstrates that the circuit is planar. Before going further, we will need some basic 
properties of planar circuits in which power is applied across a boundary edge (note: 
boundary is used here in the usual context of planar graph embeddings and should not 
be taken to mean adjacency with the sink node). 



Lemma 7.1. Let S be a planar resistive circuit with a given embedding such that the 
power source is attached across a pair of vertices on the boundary, say t (positive) and 
s (negative). There are exactly two paths from t to s along the boundary. Then the 
direction of current in any boundary edge is along these paths from t to s. 

Proof : Assume there exists an edge (a, &), in which the current flows from a to b, 
i.e. the direction opposite to the path from s to t. Since net inflow of current occurs only 
at node s, starting from a one can construct a path to node s such that in each edge 
current flows from s to a. Similarly one can construct a path from b to t where in every 
edge current flows from b to t. Because of planarity, these paths should intersect at some 
node, say c. Then we have a loop (as shown in figure (jH])) from a to 6 to c to a in which 
current flows in the same orientation in every edge. This contradicts Kirchoff's theorem 




Figure 6: Direction of current in the boundary 
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about sums of potential differences along any loop in a circuit. Hence such an edge (a, b) 
cannot exist. ■ 

This gives us the direction of current in every boundary edge, given a planar circuit 
(with some embedding) and a power source applied across a boundary edge. We now 
consider the effect of changing the resistance of some boundary edge on the potentials 
that appear on the nodes along the boundary. The following lemma tells us the change 
required in the resistance of a boundary edge to induce the desired effect on the potentials 
elsewhere on boundary. 

Lemma 7.2. Let S be a planar resistive circuit with a given embedding such that the unit 
potential is applied across an edge on the boundary, say node s (positive) and the sink t 
is maintained at zero potential. Given any boundary edge e, increasing its resistance will 
decrease the potential that appears on any node along the portion of boundary between e 
and t and increase the potentials for the portion lying between s and e. 

Proof : Using lemma mi we know that in any edge, say e, along the boundary, current 
flows in the direction of the boundary path from s to t. If we increase the resistance of 
e, i.e. dZe > (where Ze is the resistance of e), the effect on current flowing through 
any other boundary edge can be predicted using the compensation theorem. Previously 
current Ig was flowing in direction s to t and dZ,, is positive. The power source of —ledZg 
when inserted in e, induces a current in the direction i to s (again using lemma FZ.ip . 
Hence, the effect of increasing the resistance of e is that current in every boundary edge 
decrease. Since potential of any boundary node, say v, between e and t is simply ^ Zgilg' , 
where the sum is over all edges lying between t and e. Since the resistances are constant 
and currents are decreasing, the sum also goes down. The case of vertex v lying between s 
and e is analogous. Except for the fact that the potential of v in this case is 1 — ^ Ze'Ie' ■ 
Increasing the resistance Z^ decreases the summation (like in previous case) and so the 
net value increases. This completes the proof of lemma. ■ 

Continuing our discussion of the circuit S'm, we increase the resistances connecting 
any two nodes adjacent to sink to the known upper bound of |£^(>S')| — k. Using the 
lemma [7721 we know that each of these increments decreases the value of 7r„. (uj). Denote 
X = I £'(5') I — k. Figure ^ shows the circuit we have in the end. The value of TTy.{vj) 
obtained in this circuit will serve as a valid lower bound on the value we are seeking. 




Figure 7: The line circuit 



Lemma 7.3. Given the circuit Sm o,s described above. For any pair of nodes v and w 
which are adjacent to the sink, the following upper bounds on the value of TTy{w)^^ always 
hold. 

max{TT^{vy^} = Oii\E{S)\ -k + 2)'=-^) (31) 

v,w 

Proof: To keep the notation clean, we relabel the nodes in our circuit as follows. 
The node Vi is ui and Vj is Uk- AH the nodes lying in between are indexed in order of 
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occurrence on the path from Vi to vj. We apply unit potential at the site ui such that 
a potential of T:{ui) appears at node Ui, in particular 7r(wi) = 1. Next, we scale the 
potential applied at ui so that unit potential appears at node Uk- Denoting the potential 
at node Ui by Vi, Kirchhoff's equations of current conservation at any node Ui is, 

V, 
Vk 

which rearranges to give the 

V = {x + 2)V+i - V+2 (32) 
with the boundary condition 

Vk - 1 
Vk-i = x + 1 

Consider the system, V^' = {x + 2)1/'^;^, V"^! = 1. Then for each i, V{iy > V{i). Then 
V{ = (x + 2)^^'^.{x + 1). Therefore we have, 

max{^„(z;)-i} = 0{{\E{S)\ -k + 2)'=-i) 

V,W 

Note that the above solution is not far from the solution of the original set of equations. 



x + 2 



x + 1 

recursive formulation, 







' X + 2 -1 ■ 


k 


. ^2 . 




1 





The asymptotic eigenvalues of the matrix are x + 2 and (for large x) and so the 
value of Vi would be a linear combination of {x + 2)'^~^.(a; + 1) and some constant, which 
is asymptotically the same as our approximate solution. ■ 

We have already seen that. 



td{S) = 0(|£;(S')|.max{7r^(w)-i}) 

v.w 

So, for the case of a sandpile with k connections to the sink, we have the Theorem 

rm 

Remark: We observed earlier in the introductory section that the line sandpiles have 
exponential transience classes. With slight amendment, the arguments used in proving 
the bounds stated above can be used to derive exponential lower bounds on the potential 
response in the line circuit. All one needs to do is replace the value of x by 2 and reduce 
the resistance of each connection to sink to half units. A completely combinatorial proof 
of the exponential nature of the transience class of line sandpiles appears in 

8 Equivalence of Triangular and Hexagonal Sandpile 

The definition of transience class describes it as the exact number of particles which 
surely induce a toppling everywhere in sandpile. In analogy with the question of time (or 
space) complexity of algorithms which asks for the maximum time taken by an algorithm, 
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classifications exists on connected sets in C" according to the maximum possible growth 
rates of continuous harmonic functions (the classical harnack's constant) in terms of 
dimension and size of the set, upon graphs with respect to conductances, upon the speed 
of rumour spreading in graphs in terms of graph conductances [lOj , [9j , upon graphs with 
respect to the growth rates harmonic functions itself (the harnack's constant in discrete 
setting), etc. Our goal is to impose a similar classification on sandpile families. In this 
section we will show that polynomial bounds on the transience class of one sandpile can 
be used to imply polynomial bounds on a related sandpile by considering the example of 
sandpiles based on honeycomb and triangular lattices. 

An indexed family of sandpiles {S'„} is said to belong to the transience class TCL{f{n)) 
iff for all values of n 

td{Sn) = 0{f{n)) 

The transience classes TCL{exp{n)) and TC L{poly{n)) are defined in the usual man- 
ner. Our result on grid sandpiles establishes that tcl{xn) belongs to TCL{n'^). We now 
introduce the notion of transience class equivalence. 

Definition 10. We write {An} ^tci {Bn} if for any transience class TCL{f{n)), {An} G 
TCL(f{n)) ^ {B„} e TCL{f{n)). 

Two sandpile families {A„} and {Bn} are transience class equivalent if they belong to 
the same transience classes. This formalises our intent to classify sandpiles into classes, 
where the number of particles needed for complete percolation is asymptotically equal, 
upto constant factors, for every sequence. This notion assumes importance in cases, when 
a sandpile-graph can be replaced by another sandpile-graph, equivalent in the above sense 
where transience class computations are easier to deal with. We will now show that a 
family of finite sandpiles based on honeycomb lattice, say {i?„} belongs to TCL{poly{n)) 
iff the analogous family of finite sanpiles based on triangular lattice, say {Tn} belongs to 
TCL{poly{n)). 




Figure 8: Finite sections of honeycomb and triangular lattices 

In the simpler case of infinite (boundary- less) lattices, because of unbounded extension, 
it does not make sense to talk about the transience class. However, the sandpile impedence 
between any two sites is still well defined and is the right property to discuss. For planar 
lattices based on regular tessaletions of plane, these values can be estimated using simple 
particles conservation based combinatorial arguments. Given any planar lattice It is not 
too difficult to show that for any pair of vertices at a distance n (shortest path length in 
the underlying graph), the value of Rs{n) is 0{n^). The only property one needs is that 
the number of vertices in any region go up as the square of the radius of the region and 
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some symmetry properties which are integral to regular tesselations. We now consider 
the case of finite honeycomb and triangular lattices with boundaries. 

Figure [5] depicts finite sections of these lattices. The boundary edges, are connected to 
the sink node s in both cases. Consider a sequence We will construct the analogous 

sequence {T„} whose membership in TCL(poly{n)) will imply membership of {H„} as 
well. Let Hi be any member. Consider the resistive circuit based on it, also referred to as 
Hi. This circuit will be transformed into an equivalent circuit T^. In the present context, 
equivalance will have a slightly more general meaning then in electric network theory. 

Definition 11. Two sandpile circuits 5*1 and 52 with the same boundary set B{— {v\v ^ 
s}) are said to be equivalent, if for any vertex v £ B, when unit potential is applied across 
V and s, the potentials induced at all other vertices is identical in both cases. We denote 
network equivalence by Si 82- 

Following from lemmas 13.61 and 13.71 one needs to consider only the vertices in the 
boundary set for obtaining bounds on tel. Hence, when we say that the circuits Hi is 
equivalent to Ti, the bounds on tcl are identical. Since both the particles addition and last 
toppling nodes are on boundary. Lemmas 13.41 [331 and the bounds on r(.) ensure that if 
the tcl is polynomial, the bounds obtained using Theorem 13 . 1 1 are also polynomial. So for 
sandpile sequences, the polynomial transience class is closed under equivalant reductions. 

Lemma 8.1. Given {A„} and {Bn}, if Ai Bi for all values of i, then {A„} € 
TCL{poly{n)) <^ {B„} g TCL{poly(n)). 

Before we start the reducing Hi, we will need the following result. 




Figure 9: The star-delta transformation 

Property 8.1. (Star-Delta Transformation, 181) The configurations shown in figure \^ 
are equivalent for 

ab + be + ca 



and likewise for the values of B and C. 

In the context of sandpiles, as long as the central node in the star configuration is not 
in the critical set (the concerned dense subset of boundary set), one can replace the config- 
uration with the equivalent delta configuration without changing the potentials appearing 
on any boundary node when unit potential is applied at any boundary node. Consider 
the figure (TUj which demonstrates a honeycomb lattice and its equivalent triangular lattice 
superimposed in dotted lines. The star configurations belong to the honeycomb lattice 
and are made up of unit resistances. The delta configurations (in dotted lines) constitute 
the delta configuration and each resistance has value 2 units. Note that even if there 
exist boudary edges that do not belong to any complete star, they don't pose any essen- 
tial problem as every unit resistance can be replaced with two 2 unit resitors in parallel, 
as shown in the figure. The reduced triangular lattice we obtain is made up of 2 unit 
resistors. Halving each resistor's value induces a constant factor change in the values of 
7r(.) and r(.) functions over this circuit. The sandpile corresponding to this circuit is also 
denoted by Ti. We thus have a pair of sandpiles Ti and Hi such that membership of one 
in TCL{poly{n)) is equivalent to the membership of other. 
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Figure 10: A honeycomb based grid and its equivalent triangular lattice grid 



The reductions we display above, prove important in the cases when only one of the 
members of an equivalent pair has the necessary symmetries to deduce polynomial bounds. 



9 Future work and Open problems 

The main open question is that of tightening the bounds on te^(GRID„). As noted in the 
remark at the end of subsection 15.51 one can expect substantial improvements only in 
the estimation of max^, Trw{v)~^. We believe the approach using Lemma |6. II is the most 
promising avenue. 

The general form of the eigenvalues and eigenvectors of a grid are well known. Using 
these and the results in Lemma 16.11 one can approximate (up to constant factors) the 
corner to corner potential correlation using the following function. 



n\ ^ 4_2cos(^)-2cos(^) 



After spending considerable time on trying to resolve the absolute size of this expres- 
sion, which started as a seemingly harmless looking question and subsequently led to a 
formulation of sorts on the general size estimation problem of alternating sums based on 
uniformly continuous functions, the authors must admit their inability in resolving this 
rather technical problem and invite the interested reader from the theory community to 
take it up from this expression. For purposes of restricting the manuscript size, we have 
not described the complete algebra leading to this expression as well as the details of our 
estimation procedure, these are available on personal request. 

The second open question is showing that the opposite corner sites are indeed the worst 
pair for single site particle addition strategies. The more general question of finding the 
analogous pair in general graphs is also interesting. We conjecture the following implicit 
characterization of such a pair. 

Conjecture 2. For a given sandpile S , if the site u allows one to attain the worst case 
bounds for single site particle addition strategies, then in the corresponding circuit, there 
exists a boundary site w such that tTu]{v) — minp^qiTp^q) . 

This hints at another of the many ways in which sandpile are similar to electric 
networks. We have already demonstrated that single site particle addition strategies are 
not enough to attain the transience class. We believe however that the following weaker 
conjecture holds. 

Conjecture 3. Assume that for sandpile x, the induced sub-graph on the set of ordinary 
vertices is connected. Then the transience class of x ( the largest weight of any transient 
configuration ) is bounded from above by the sum of height of the tallest transient stack of 
grains placed on a single site and the size of graph x- 
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Remark: One can show that the transience class of a sandpile is bounded from below, up 
to constant factor, by the graph size. This follows from the fact that when every site has 
toppled, no two adjacent sites can both have zero particles. Which implies that for every 
edge, at least one particle is on board. So tcl{x) = ^{\x\)- In the light of this observation, 
conjecture [3] means that the estimates derived using single site particle addition strategies 
are constant additive factor approximations of the actual transience class. Also note that 
the above conjecture is stronger then the question raised by Sunic, mentioned in [31, that 
local (i.e. single site particles addition type) transience classes are bounded by 0.5-factor 
approximations of the transience class, as we have an additive error term compared to 
the previous multiplicative one. 
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